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Abstract: We explore entanglement entropy of a cap-like region for a generic quantum 
field theory residing in the Bunch-Davies vacuum on de Sitter space. Entanglement entropy 
in our setup is identical with the thermal entropy in the static patch of de Sitter, and we 
derive a simple relation between the vacuum expectation value of the energy-momentum 
tensor trace and the RG flow of entanglement entropy. In particular, renormalization of 
the bare couplings and logarithmic divergence of the entanglement entropy are interrelated 
in our setup. We confirm our findings by recovering known universal contributions for a 
free field theory deformed by a mass operator as well as obtain correct universal behaviour 
at the fixed points. Simple examples of entanglement entropy flows are elaborated in 
d = 2,3,4. In three dimensions we find that while the renormalized entanglement entropy 
is stationary at the fixed points, it is not monotonic. We provide a computational evidence 
that the universal ‘area law’ for a conformally coupled scalar is different from the known 
result in the literature, and argue that this difference survives in the limit of flat space. 
Finally, we carry out the spectral decomposition of entanglement entropy flow and discuss 
its application to the F-theorem. 
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1 Introduction 

Entanglement entropy is an important tool in condensed matter physics [1-4], quantum 
held theory [5-9] and quantum gravity [10-21], Nowadays this technique is in the spot 
light of nonperturbative studies of the structure of quantum held theories (QFT). In 
particular, constraints imposed by c-theorems, which describe the irreversible nature of 
renormalization group (RG) hows between the hxed points, is one of the examples where 
entanglement entropy methods are applied [22-33]. 1 These theorems spring directly from 
Zamolodchikov’s proof [42] of irreversibility of RG hows in two dimensions. In general di¬ 
mension, however, the mechanism of irreversibility is still unclear, but Cardy’s conjecture 
[43] suggests that this feature is inherent to any QFT for any even dimensional space- 
time. Recently, a striking proof was found in four dimensions [44], see also [45, 46], but its 
generalization to six and higher dimensions is not straightforward [47]. 

1 See also [34-41] for recent progress in applying entanglement entropy techniques in a QFT context. 
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A natural extension of Cardy’s conjecture for quantum field theories in an odd di¬ 
mensional space-time was proposed in [23, 24], Based on the holographic studies of en¬ 
tanglement entropy in general dimension, they suggested that the universal coefficient of 
entanglement entropy for a spherical region is decreasing along the RG trajectory. At the 
fixed points of RG flows in even dimensions this coefficient is proportional to the central 
charge used by Cardy to formulate his conjecture [6, 48-50] , while for odd dimensions it 
was shown in [50], see also [51], that this coefficient is directly related to the T-theorem 
[52, 53]. These observations indicate that entanglement entropy provides a useful frame¬ 
work to study c-theorems in general odd and even dimensions. 

Important progress in understanding c-theorems in three dimensions was made in [54] 
where Casini and Huerta used strong subadditivity of entanglement entropy [55] to prove 
the three-dimensional F-theorem for any unitary and Lorentz invariant QFT. There is not 
yet an alternative derivation. However, it is certainly important to understand the key 
insights of the proof in terms of conventional held theoretic methods, and one of our goals 
in this paper is to make a step in this direction. 

Motivated by [56], we think of entanglement entropy as a scalar functional of the held 
theory couplings and geometry of the setup. Imposing various symmetries inherent to the 
underlying QFT yields a set of Ward identities and RG how equations satished by the 
entanglement entropy functional [57]. The resulting identities depend on various correla¬ 
tors involving insertions of the modular Hamiltonian, and therefore from computational 
perspective they are not really tractable since the modular Hamiltonian is not known and 
not local in general. However, for a special class of co-dimension two entangling surfaces 
which exhibit rotational symmetry in the transverse space, the modular Hamiltonian can 
be identihed with the angular evolution operator around a given entangling surface [58-60] 
while the corresponding Ward identities and RG how equations can be written in terms of 
the standard correlation functions. 

Flat entangling plane in Minkowski space-time is probably the simplest realization of 
such a symmetric setup. It has been shown in [60] that entanglement entropy in this case 
equals thermal entropy of the Rindler wedge. Perturbative studies of this setup are useful 
for understanding the structure of both the universal entanglement entropy [56, 61, 62], see 
also [63-65], and the modular Hamiltonian [66]. There is also a way to use it in order to set 
a precise relation between the entanglement entropy and renormalization of the Newton’s 
constant [33]. However, this setup has no built-in scale, and its topology is trivial, which 
makes it hard to believe that one can learn much about the non-perturbative aspects of 
c-theorems in higher dimensions based on this example. 2 

Our proposal in section 2 is to use a rotationally symmetric entangling surface for a 
field theory residing in the Bunch-Davies vacuum on de Sitter space. 3 4 Rotational symmetry 
ensures simple structure for the modular Hamiltonian while the curvature of maximally 
symmetric geometry effectively sets the renormalization group scale. 1 Such a setup is 

^However, see [33] for an alternative derivation of Zamolodchikov’s theorem using this setup. 

3 For a CFT this construction was studied in [24] and holographically in [23, 24], see also [67] and [68] 
for various aspects of the entanglement entropy in de Sitter space. 

4 In the system we study, the scale of entangling surface is also defined by the radius of de Sitter space, 
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scalable in the sense defined by [26], and RG flow equation for entanglement entropy can 
be expressed in terms of two-point function of the energy-momentum tensor. A further 
simplification takes place after we notice in section 3 that the RG flow of entanglement 
entropy has simple interrelation with the vacuum expectation value (vev) of the energy- 
momentum tensor trace. This leads us to a quite interesting conclusion that the logarithmic 
divergences associated with renormalization of the bare operators and bare parameters 
in the energy-momentum tensor on de Sitter space are simply related to the universal 
entanglement entropy. We illustrate this peculiarity in section 4 using the example of 
massive free field theories, where renormalization of the cosmological constant must be 
done even in the absence of interactions. 

Our construction has much in common with the Rindler space and with the observa¬ 
tions made by [60]. For example, entanglement entropy equals thermal entropy for any 
QFT in our setup. However, there are several distinctions which make it more interesting. 
Firstly, the Euclidean continuation of the causal domain of entangling region in de Sitter 
space has non trivial topology in comparison to its counterpart in Minkowski space (sphere 
versus flat plane). Secondly, the curvature of de Sitter space sets up a characteristic tem¬ 
perature which manifests itself through the finite terms in entanglement entropy. Finally, 
one can explore an RG flow of entanglement entropy between the fixed points by varying 
the radius of de Sitter space. Given that de Sitter space is more realistic than Minkowski 
background, such RG flows are worth investigation. 

In section 4 we build on the advantages of our setup to study simple RG flows. In par¬ 
ticular, we find that renormalized entanglement entropy defined in [26, 54] is not monotonic 
on a three dimensional de Sitter space. Of course, this observation does not contradict [54] 
since there is no reason to expect that RG flows on de Sitter and Minkowski space-times 
share the same merits. Furthermore, in contrast to the flat space results in [69], our findings 
in section 4 indicate that the renormalized entanglement entropy is stationary at the fixed 
points on a sphere. This qualitative difference between the backgrounds can be attributed 
to the absence of infrared (IR) divergences in our setup, since sphere can be thought of as 
a finite box. It should be noted, however, that this result is in certain tension with [69], see 
also [70], where non-stationarity was found not only for a scalar field theory living in flat 
space, but also for a conformally coupled scalar on a sphere. This disagreement is closely 
related to another observation made in section 4 where we argue that the universal ‘area 
law’ for a massive conformally coupled scalar field on a sphere is different from its counter¬ 
part in flat space [71-73]. This discrepancy is not sensitive to the size of de Sitter radius 
which points in favour of controversial [15, 33, 36, 74-78] difference between the universal 
entanglement entropies for minimally and conformally coupled scalar fields in flat space. 
Possible explanation of this phenomenon is given in [79], see also [56, 62, 80]. 

Finally, in section 5 we use [81] to elaborate the spectral decomposition of entanglement 
entropy flow in general dimension. This decomposition is determined by an integral over 
a positive weight function which represents the spin-0 part of the spectral representation 
for two point function of the energy-momentum tensor. In particular, the rate of change of 

therefore there is only one global geometric scale. 
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entanglement entropy with respect to RG scale has definite sign, provided that the integral 
over the weight function converges. In three dimensions it suggests that the finite part of 
entanglement entropy may play a role of the c-function. However, this conclusion is too 
naive. Indeed, possible logarithmic divergences of the entanglement entropy for a QFT 
in three dimensions [56, 82] is a clear sign of potential problems with convergence of the 
spectral integral. Moreover, finite part of entanglement entropy at the IR fixed point is 
contaminated by various mass scales which are remnants of RG running into the fixed point 
[33]. In the deep infrared these finite terms are part of UV physics, and therefore they 
should be subtracted to distill the genuine universal piece of entanglement entropy at the 
IR fixed point. 


2 Setup 


Let us consider a generic QFT living on a d-dimensional background de Sitter space of ra¬ 
dius R- a submanifold of (d+l)-dimensional Minkowski space described by the hyperboloid 
of one sheet, 

x\ + X 2 + ••• + xji — Xq = R? . (2-1) 

We will assume that the field theory resides in a unique Bunch-Davies vacuum state, also 
known as the Euclidean vacuum, which is invariant under all the isometries. In global 
coordinates, 


xq = R sinh ( — 
- R 


Xi = Rcosh j, i = 1,.., d , 
with L2? = 1, the metric on de Sitter space takes the form 


ds = — dt + R cosh ( — )dQ d -i 


( 2 . 2 ) 


(2.3) 


In what follows we aim to consider entanglement entropy of two equal cap-like regions 
A and B which are located on a xq = 0 (f = 0) time slice of dS d . This slice is invariant 
under t —> — t which makes it possible to analytically continue the metric to the Euclidean 
signature, i.e., xq — > ixd+i or equivalently t —> it. In particular, the problem of finding 
entanglement entropy boils down to studies of a QFT living on a d-dimensional sphere of 
radius R, see Fig.l. 

By construction, the entangling surface, E, divides the equator of S d into two equivalent 
cap-like regions. In particular, it exhibits rotational symmetry in the transverse space 
which has topology of a 2-sphere. Hence, we find it convenient to foliate S d such that this 
symmetry becomes manifest, 


Xi = R sin 6 for 1 < i < d — 1 ,O<0<^- 

Xd = R cos 9 cost , 

Xd +1 = R cos 9 sinr , 




1=1 


(2.4) 
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Figure 1. Illustrating the setup of a QFT living on an S d sphere. The entangling surface (at t = 0 
and 9 = 7r/2) divides the equator into two equal cap regions (shown in blue and red). 


where (r, 6) parametrize a two-dimensional transverse space to the entangling surface at 
9 = j. 5 The infinitesimal line element on S d takes the form 

d +1 

ds 2 = ^ dx 2 = R 2 ^ cos 2 Odr 2 + dO 2 + sin 2 9d^L 2 d _^j . (2.5) 

2—1 

It worth noting that the above foliation of S d corresponds to analytic continuation of 
the static patch of de Sitter to Euclidean signature, and it has been shown in [50] that 
if the field theory is conformal, then entanglement entropy for a sphere of radius R in 
Minkowski space is equivalent to the thermodynamic entropy of the thermal state in the 
static patch of de Sitter geometry. Of course, in the absence of conformal symmetry there 
is no simple relation between these entropies. However, by construction it is still true 
that the entanglement entropy of two equal cap-like regions for a theory residing in the 
Bunch-Davies vacuum equals to the thermal entropy associated with a thermal state in the 
static patch of de Sitter. A quantitative manifestation of this (exact) relation is provided 
in section 4. 

Now let us assume that the field theory on S d is some CFT deformed by a set of relevant 
operators Oi of scaling dimension A* and associated couplings A*. Since entangling surface 

5 In two dimensions — ^ < 9 < and entangling surface consists of two disjoint points at 9 — 
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exhibits rotational symmetry in the transverse space, the modular Hamiltonian is given by 
K = -2vr [ T^en v + d , ~ ' ' 2 51 


y/gSg^ix) ’ 


( 2 . 6 ) 


where / is the Euclidean action of the field theory, the integral runs over the region A 
of the equator, n v = (i?cos0) _1 9 r is normal to this region, = d r is the Killing vector 
associated with rotational symmetry around £ and d is an additive constant to ensure 
proper normalization of the reduced density matrix, i.e., Tre -A = 1. 

Now recall the general flow equations that describe changes in the entanglement en¬ 
tropy under flow of the relevant couplings and deformation of the geometry [57], 
dS EE f . 1 61 


dX 

SS EE 


f = -/Av7W(0.(x) K) X , O(x) = 
\/g(x) 


-{T^x,)K) X 


(2.7) 


( 2 . 8 ) 


5g^ u (x) 2 

where A collectively denotes all the couplings. These flows are directly related to the 
first law of entanglement [83-86] 6 . Of course, for a general deformation K will not be 
given by (2.6) since rotational symmetry will be destroyed. However, if we restrict our 
consideration to a particular one parameter family of constant rescaling of the background, 
then rotational symmetry is preserved and (2.8) takes the following form 

/ Ar s^%) = J iix OMCr(x)A-h , (2.9) 

where T( x) = g^ v T^ u {x) is the trace of the energy-momentum tensor, and integrals run 
over S d . 

The left hand side simplifies if we notice that for our choice of the foliation constant 
Weyl rescaling of the background metric is tantamount to a constant rescaling of the radius 
of S d , and we obtain 

.dS EE 


R- 


dR 


= J d d x y/g(x)(T(x)K)x ■ 


( 2 . 10 ) 


Eq. (2.10) can be also derived using the standard replica trick approach just because 
the rotational symmetry around the entangling surface makes it possible to introduce a 
well-defined conical defect, such that 


S EE = lim I + 1 ) log Zi_ e , 
e-s-o V de 


( 2 . 11 ) 


where Z\_ e is the partition function on the n-fold cover of a sphere and e = 1 — n. Hence 
[13, 24], 


R- 


. dS,, 


dR 


d 

= -2 lim I — + 1 
£—V de 


d 

= lim I — + 1 
e—V de 


d d xg^(x)- 


6g^ v (x 
d d x y/g(x){T(x)) i_ e 


■ log Zi_ e 


= / d d x yjg(x){T{x)K) 


( 2 . 12 ) 


6 See also [87] for a discussion of certain parallelism between the entanglement and the laws of thermo¬ 
dynamics. 
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where the last equality rests on [37, 88] . 

Unfortunately, (2.10) is ambiguous when supports of the energy-momentum tensor and 
the modular Hamiltonian collide. Hence, we resort to the Ward identities which help to 
circumvent this ambiguity. 

3 Ward identities 

In this section we use the Ward and trace identities to constrain correlation functions 
involving the modular Hamiltonian and the energy-momentum tensor. They will help us 
to fix certain ambiguities involving situations when the energy-momentum tensor collides 
with the support of modular Hamiltonian. 

Let us consider the vacuum entanglement entropy of a generic surface E for a quantum 
field theory defined on an arbitrary curved space. Assuming that regularization preserves 
invariance under the diffeomorphisms g pL , —> where and v p stand for a 

covariant derivative and arbitrary vector field respectively, and taking into account that 
entanglement entropy is a scalar functional, implies' 

o = ~ J d d x{V»v v + VV*)^J^y^ EE = - J d d x sfg v u \7 fl (T flu (x)K) , (3.3) 

where in the second equality we substituted the definition S EE = —Trplog/5 = (K) and 
used normalization of the density matrix to drop [5K/5g^ u ) = 0, see [56]. Hence, we 
conclude 

V»(T^{x)K) = 0 . (3.4) 

Let us contrast this formula with the standard Ward identity for the two point function 
of the energy-momentum tensor. Starting from the counterpart of (3.4), 

V^V(x)) = 0 , (3.5) 

which one obtains from the condition that the effective action, W = log Z with Z being a 
partition function, is a scalar functional under diffeomorphisms, and differentiating it with 
respect to the metric, yields [81] 

V»(T^(x)T a p(y)) = Vj, (S^6 d (x,y)) {T ap (x)) + 2V CT (s^ 6 d (x, y)(T pu (x))) , (3.6) 

where 5 d (x,y) = 5 d (x — y)/ y g(x), whereas the 2-point function is defined by 

. . 5 2 W 

(T^(x)T a p(y)) = 4 6g ^ x)Sga p^ • ( 3 - 7 ) 

'The most general Ward identity of this type, which also accounts for the possibility of local sources 
takes the form 

J d * x (-^v I '+V , 'v») I ^+v»d,X^)S EE = 0, (3.1) 

or equivalently 

V^(T^(x)K) = d^iOi^K) . (3.2) 

Ultimately, we will be interested in A ! ( x ) = A 1 the physical coupling constants, and therefore we drop terms 
involving 9 M A\ 
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The right hand side of (3.6) vanishes up to contact terms involving ^-functions, while 
the right hand side of (3.4) vanishes identically. In particular, in the special case when 
the modular Hamiltonian is local and is given by (2.6), it would be naive to use (3.7) to 
evaluate the right hand side of (2.10). One needs to modify this correlator in the limit of 
coincident points such that (3.4) holds identically. A necessary modification in a slightly 
different context was carried out in [81], and we review it here. 

One starts from noting that the vacuum expectation value of any local scalar operator 
on a sphere is just a constant and also 

( T /lu (x )) = — — g^vix) , ( 3 - 8 ) 

where C is some dimensionless function of the couplings. Eq. (3.5) is trivially satisfied, 
whereas (3.6) becomes 

V M (T^(x)T 0 , / 3 (y)) = ~2Jtd { s/ ^ d (x,y)g af) {y) + d d (x,y)g pu (x))^j (3.9) 

Next define, 

1 C 

{Tf_n/(x)T a /3(y))con = {T)u/(x)T a p(jj')) + — 4“ 9fi/39i'a T 9^iv9afi)d {%->y) ■ (3.10) 

where ( ) con means “conserved”, as in [81]. This 2-point correlator satisfies a desired 
conservation equation [81], 

V^(T^(x)T aP (y)) 

con — 0 i (3.11) 

and therefore it should be used in (2.10) to ensure that the Ward identity (3.4) is satisfied. 
Now recall that for our setup homogeneous Weyl rescaling obeys 

In particular, applying this operator to the both sides of (3.8) leads to 
- J d d y^/g{y)(T(y)T^(x)) + d(T pu (x)) = * r£ r g^{x) - ^ Rd g^(x) . (3.13) 

Or equivalently, using (3.10) 

- J d d yV9¥){T(y)T^(x)} m = R jL (f7) , (3.14) 


Substituting this formula into (2.10) and using (2.6), we obtain one of our main results 8 


R- 




2ir „ d ( C 


dR 


d R dR\R d ) JJ' n 


n H R d+1 d ( C 


d dR V R d 


(3.15) 


8 We stress that unlike the one point function ( K }, the (divergent) constant d drops out of connected 
correlators. 
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where is the surface area of a unit d -dimensional sphere, 



It worth mentioning that (3.15) relates finite quantities on both sides of the equation 


since C is finite by definition. In particular, this formula eliminates the power law diver¬ 


gent terms of EE such as the well known ‘area law’. Of course, there is nothing bad about 
absence of these terms since they are scheme dependent and therefore vanish under appro¬ 


priate choice of the regularization scheme ( e.g dimensional regularization). However, the 
logarithmic divergence of EE is always retained, and its coefficient is scheme independent. 
This observation may cast doubts on possible applications of (3.15) for a QFT in even di¬ 
mensions, where the universal EE is associated with a logarithmic divergence. In section 4 
we elaborate on simple examples which clarify this subtlety. We show that the universal EE 
in our setup is directly related to the logarithmic divergences of the bare parameters defin¬ 
ing the theory, and the standard renormalization results in a finite EE with the universal 
data being encoded in the coefficient of the (finite) logarithmic running. Throughout the 
paper we apply dimensional regularization scheme, and therefore non-universal divergences 
of EE will be absent from our calculations. 

3.1 Thermal interpretation 

As pointed out in section 2, entanglement entropy equals thermal entropy in the system 
under study, S EE = S Th . Therefore (3.15) should have a simple interpretation in terms of 
standard thermodynamics. We provide such an argument below, see also [50] for thermal 
analysis of S EE = S Th at the fixed point. 

Let us start from the following thermodynamic relation, 


S Th = P(U - F) 


(3.17) 


where U and F are thermal and free energies respectively, and (3 = 2 irR is the inverse 
temperature. As usual, U is given by the expectation value of the operator which generates 
translations around the thermal circle parametrized by r in (2.5). Hence, based on (3.8) 



(3.18) 


Now we can operate with R4^ on both sides of (3.17). Using (3.12), our convention (2.6) 
for the energy-momentum tensor, and recalling that / 3F = —W, yields 



(3.19) 


This is exactly (3.15), thus we showed that S' EE = S' Th . Since at the fixed point we have 
See = SThi we conlcude that this is true also outside the fixed point. 
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3.2 RG equation for entanglement entropy 


The renormalization group running of entanglement entropy in our setup can be readily 
evaluated. Based on (3.8), we have 

% = - Rd ^ {T{x)) = Rd J dd yV^y)my)nx)) + dR d (o t (x)) = (*)>) , 

(3.20) 

where in the last equality we used homogeneity of the sphere to replace integral over y with 
integral over x and then applied the correspondence (3.12). Next we combine this result 
with the trace Ward identity 9 [81, 89] 

(T) + (d - A - P)\(0) = A => C = R d ({d-A-/3)\(0) - A) , (3.21) 

where /3 is the beta function, and A is the trace anomaly formed from the Riemann ten¬ 
sor and its derivatives, i.e., A = 0 for any odd dimensional space-time while for even 
dimensional sphere A oc R~ d . As a result, we obtain [81] 

(r±-V-A- «4) c - -R± (/fU) - 0 . (3.22) 

This equation combined with (2.7), (3.8), (3.15) and the definition of the modular Hamil¬ 
tonian (2.6) leads to A|^ = — ^f^A^ + &dR d (O(x)), 10 substituting this we get: 

( U^-{d-A-p)\-^jS EE = n d R d A . (3.23) 

Given that A is constant on a sphere, we deduce that the right hand side is just the 
integrated trace anomaly. This equation is a particular realization of the general idea 
presented in [57]. Here, however, we keep the finite anomalous term on the right hand 
side of (3.23) since we are interested to evaluate finite logarithms which combine with the 
universal divergence of entanglement entropy to form a dimensionless term. 

3.3 Conformal fixed point 

Let us consider a rotationally symmetric entangling surface for a CFT residing in a vacuum 
state on some curved manifold. 11 In this case, the right hand side of (2.9) is completely 
fixed by the trace anomaly 

{T{x))=A = J2 b nIn(x)-2{-l)iaE d (x) + B , V lx J^x) , (3.24) 

n 

which defines the central charges for a CFT in an even number of dimensions. Each term 
on the right-hand side is constructed from the background geometry. I n are Weyl invariant 
combinations of the Weyl tensor, the Cotton tensor and the Bach tensor as well as their 

9 For brevity we suppress index i of the relevant couplings. Sum over this index is implicitly assumed. 

10 When deriving this relation it should be noted that the difference between A^ and A stems from 
the fact that ^ ~ g^ (O) but (|jy) = 0, as can be seen from the requirement Tre _if = 1 

J1 It should not necessarily be our setup, e.g., waveguide geometry [71, 73] works equally well. 
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covariant derivatives. These basis tensors all vanish on any conformally flat background. 
The last term in eq. (3.24) is a scheme-dependent total derivative which can be eliminated 
by adding a covariant counter-term to the effective action. Finally, E^ is the Euler density 
in d dimensions, 


E2p(E) 


1 

(87r)p r(p + 1 ) 


sy 1 V2p—1 V2p 2 

Ufll 112 M2p-1 /V ^1^2 


'JP ^2p-1^2p 
/V ^2p-1^2p 5 


(3.25) 


where 6^ )) 2 2 ■■■denotes a totally antisymmetric product of 2 p Kronecker delta sym¬ 
bols and the normalization ensures that f sd d d Xy/g Ed = 2. 

Varying (3.24) with respect to g tll '{y) and using the definition (3.7), yields 


( T (x)Tpu(y)) ~ 6 d (x,y)g flu (x)(T(x)) - 2 8 d (x,y){T flv {x)) = 


2 5A{x) 

Vgiyjta^iy) ' 


(3.26) 


In our setup the vev of the energy-momentum tensor takes a simple form (3.8), and therefore 
according to the definition of conserved 2-point function (3.10), we may identify the left 
hand side with {T(x)T fll ,(y)) con . We will assume that similar redefinition exists for any 
system with rotationally invariant entangling surface. 12 Hence, (2.9) at the fixed point 
reads 

= - 2 "/*v 5 w/Av 51 ■ < 3 - 27 > 


where denotes the induced metric on a region A enclosed by the entangling surface. 

The right hand side can be readily evaluated in our case. Indeed, I n terms in the trace 
anomaly (3.24) play no role since they are at least quadratic in the building blocks which 
vanish for a conformally flat background. Total derivatives can be ignored, 1 ' 1 whereas the 
contribution of the Euler density is easy to evaluate since its integral is a topological term, 
i.e., 


° = Sg^hy) / ddx ^^ Ed ^ = J dd ^ " ^g^(x)E d ( X )6(x - y)^j . 

(3.29) 

Hence, 

J y = N// ^ 9pv(y) E d (y) . (3.30) 

12 It would be interesting to work out the details of this redefinition based on the requirement that Ward 
identity (3.4) holds. Our assumption here rests on the observation that in the replica trick approach there 
are no (5-functions away from the entangling surface. 

13 Variation of the total derivative is given by 

<5(V M E) = (<5V m )E + ■ (3.28) 

First term on the right hand side vanishes since = 0 on a sphere, whereas the integral of the last term 
vanishes since sphere has no boundaries. Therefore we conclude that the net contribution of the total 
derivatives to entanglement entropy flow is zero. 
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Now Ed is obviously constant on a sphere, therefore we can use our choice of normalization 
condition to get 

E d \s* = • (3-31) 

Combining, yields 

a f^y^W) t-n = 4(-l )i- l a . (3.32) 

In accord with [23, 24], see also [50]. 


4 Free fields on sphere 

In this section we use free massive fields to elaborate on various properties of the formalism 
that we have developed in the previous section. 


4.1 Conformally coupled scalar 

Let us consider a free massive scalar field on a d-dimensional sphere of radius R 

1 = J sd {\^ 2 + \ m2 ^ + ’ ( 4J ) 

where £ c = ^d-i) the conformal coupling, 1Z = is the Ricci scalar of a sphere and 

m 2 is the mass of the field. The energy-momentum tensor is given by 


Tfxu — g^ w 


(^(<V 2 + \™ 2 4> 2 + ^ c R(j> 2 ^ 


+ + £e (57^V 2 


Hence, 

T = g ,lv T ilv = -mV - ^<X-V 2 + tell + m 2 )(j) . 
Discarding equation of motion operator leads to 


v^v.) 4 ? • 

(4.2) 

(4.3) 


V = m 2 l? d V) 


(4.4) 


The vev of 4> 2 is given by the coincident point limit of the Green’s function which 
solves the Green’s equation on S d 


1 


d 


R 2 sin rf 1 x &X 


sin d 1 ) + icR + m 2 G m (x) = 5 d (Rx) 


(4.5) 


where x is the polar angle on a sphere, and we used rotational symmetry to bring one of 
the points to the north pole. 

To solve the above equation, one needs to impose regularity at % = 7r, and 


GAx) ~ W for x<<1 ' 


(4.6) 
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which corresponds to a scalar potential created by a unit charge placed at x = 0- The 
general solution which satisfies these conditions is 


_ R 2 ~ d r(A)r(-A + d-i) n , d 2 Xx 

Gm ^ (4vr ) d / 2 r(d/2) 2 Fi(\ , d 1 A , 2 , cos 2 ) , (4.7) 


A = d — ] L + iJ( m R) 2 - ^ . 


(4.8) 


Taking the limit x ~^ 0 and dropping a mass independent (divergent) constant 14 , we have 

(4.9) 


(?) - sin (£(« - 2A)) . 


7r(47r) rf / 2 i? rf 2 v 2 

In particular, (cj) 2 ) is finite for odd d, whereas for even d it exhibits a simple pole which 
corresponds to the logarithmic divergence. We treat these cases separately. 


Odd d 


For odd d, we have 


& = (T; H m - 1/2-2? - 1/4 + m 2 i? 2 ) 

(47T) / 1 ( 2 J . 1 


^m 2 R 2 - 1/4 
TT(mR) 2 coth( 7 Ty / m 2 i ? 2 — 1/4) (—)~ 

yJm 2 R 2 — TJZ (4vr) rf / 2 r (0 

-l , (d-l)(d 2 -5d + 3) d _ 3 


(4.10) 


X (: mR) a ~ l + 


24 


m,R) c 


...+ 


7r 


2 r (^)r(i^) 


Appearance of the hyperbolic function in both the C$ and the entanglement entropy after 
using (3.15) is not surprising. Indeed, such functions are the direct associates of any 
thermal computation, whereas the state in de Sitter space has thermal interpretation. To 
get rid of the thermal effects, let us consider the behaviour of entanglement entropy in the 
IR limit, mR 1. In this regime curvature corrections are negligibly small, while thermal 
effects are exponentially suppressed, and we get 


d -1 

C = (~) 2 * 
^ mH» 1 (4 7r )<i/2p 


(mR) d + d[d ~ 2)( : d ~ A \ mR ) d - 2 


24 


+ ... 


J Qscalar 
Ct-O un iv 


d +1 


(d — 2 )(d — 4) (—) 2 7 r 
dR mi?»i 24(d — 1) ( 47 r )^rW^ 


A^rri 


d—2 


+ ... 


(4.11) 

(4.12) 


where ellipsis encapsulate curvature corrections to the leading order term also known as 
universal ‘area law’. As shown in Appendix A, this term is identical to the universal ‘area 
law’ of entanglement entropy of a half space for a conformally coupled scalar field residing 
in the Minkowski vacuum. This result is a consequence of the fact that any surface and 
any background are locally flat. 

14 The mass independent term behaves as x -2 ~ , an< A therefore it diverges as power law in d > 2. Such 
divergences are scheme dependent, and therefore one can choose a particular scheme where it vanishes, e.g., 
in dimensional regularization one gets X 2 ~ d —•► 0 f° r d < 2, and therefore analytic continuation to higher d 
also vanishes. 
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Renormalized Entanglement Entropy in 3D 

For a vacuum state in 3-dimensional flat space QFT, the general pattern for EE of a disk 
is 

‘S'ee = ci — — Co , (4-13) 

where R is the radius of the disk and 6 is a UV cut off. At the fixed point, co and ci 
are some constants , 15 whereas outside the fixed point Co is a function of R and various 
scales characterizing a given QFT. Obviously, ci is scheme dependent while Co is universal, 
therefore the authors of [26, 28, 54] defined the so-called renormalized entanglement entropy 
(REE) to isolate the universal contribution , 16 

7 O 

s 3 = R -^r - See = RS' ee - S EE . (4.14) 

Using strong sub-additivity of EE [55], Casini and Huerta proved that S EK < 0 for 
a disk in flat space. Hence, = RS EE < 0, which clearly indicates that REE has a 
monotonic RG flow in flat space, and [54] 

/*oo 

Ac 0 = C - c[f = - / dRR S" E > 0 . (4.15) 

Jo 

Of course, the RG flow itself is ambiguous but the fixed points satisfy an F-theorem 
[52, 53]. Furthermore, conformal symmetry relates the REE at the fixed points in flat space 
to the universal EE in our setup [50], hence it is worthwhile to explore REE on a sphere 
outside the fixed points. 

For a conformally coupled scalar held, S EE can be readily evaluated based on (3.15), 
(4.4) and (4.9). Fig. 2 shows the corresponding plot, and it can be seen that S'ee changes 
sign. A similar issue is discussed in [53] for the free energy on a sphere, where the authors 
argue that additional subtractions are necessary because of emergence of the cosmological 
constant term in the IR limit. Although REE is not monotonic, other choices of subtraction 
scheme might result in a monotonic how as was suggested for the free energy in [53]. 

Moreover, REE in our setup exhibits stationarity at the hxed points, i.e., S EE —> 0 as 
mR —> 0, oo. This result should be contrasted with numerical studies in [69], where it has 
been shown that REE for a disk in hat space is not stationary in the massless limit. Similar 
considerations for the free energy showed that although the IR divergent free energy obeys 
stationarity, the existence of a subtraction scheme that maintains stationarity, monotonicity 
and analyticity is not obvious [53]. 

We have calculated Aco ~ 0.0638 which agrees with the expected result from the free 
energy calculation of a conformally coupled massless scalar in the UV and an empty theory 
in the IR, in spite of thermal effects mentioned above. 

15 At the IR fixed point there could be additional terms, which are remnants of the RG flow [33], e.g., the 
universal ‘area law’ ~ mR for massive QFTs [71-73]. However, the characteristic (relevant) scale of such 
terms is very large in the deep IR., and therefore in this regime they are not really distinguishable from the 
R/5 term. Hence, we do not write them out. 

16 See [26, 28] for definition of REE in general dimension. 
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d 2 S 

d(mR) 2 



• d? S 

Figure 2. d ^ m ^2 for a massive conformally coupled scalar field on a 3-sphere. Unlike its counter- 
part in Minkowski space, changes sign around mR ~ 1.6. 


Even d 

For even do the pole structure of C$ is given by 

r , yt ,/2 2„2 n$ir 2>/2 (( rf o /2 - 1/2 - i ) 2 + in 2 fl 2 - 1/4) 2 

W-l J " .... ^ d-d 0 


+ 


(4 7r )d 0 / 2 r 

(47r) d °/ 2 r J V 24 

d 0 {d 0 - 2)(d 0 - 4)(d 0 - 6 ) (5d§ - 18d 0 + 4) , _ 

- (mR) u 


+ ... 


5760 


+ ... + 2 


r 2 (f) 

do — 2 


(mR) 4 


(4.16) 

2 

d — do 


The universal divergence of the entanglement entropy is obtained by substituting this 
expression into 1 ' (3.15) 


R 


ds°: 


(do — 2) (d 0 — 4) (_)<W2 +i Al 


dR 


12(d - 1) (47T-)( rf 0—2)/2p 


x [ m d ° 2 + 


(do — 6 ) (5dg — 18do + 4) rn d ° 4 


120 


i ? 2 


+ ... + 24 


r 2 (f) m 4 

d 0 (d 0 - 2) 2 


(4.17) 


log(md) . 


The leading order term in the limit mR 1 represents the universal ‘area law’. In 
this limit curvature effects can be disregarded, and the system can be approximated by 
a free scalar field theory living in the Rindler wedge [90], therefore it should be possible 
to recover the same result by considering the entanglement entropy of a half space for 

17 As usual, simple poles in the dimensional regularization scheme correspond to logarithmic divergences, 
and we use the following dictionary log(m5) = d _) d(j , see Appendix A. 


- 15 - 






















a massive scalar field residing in the Minkowski vacuum state [60]. We carry out this 
computation in Appendix A and find full agreement. However, the resulting expression 
for the universal ‘area law’ is not the same as in, e.g., [71-73]. Possible interpretation 
of this discrepancy is given in [79], see also [56, 62, 80]. In Appendix B we present an 
independent computation of the curvature correction to the universal ‘area law’ and find 
full agreement with (4.17). Note also that unlike odd d, the universal entanglement entropy 
in even dimensions is not affected by de Sitter temperature. Absence of thermal corrections 
is due to the fact that UV divergences of EE are state independent, see e.g., [82]. 

So far we discussed how the logarithmic divergences of EE are encoded in (3.15). In 
particular, we used a naive bare expression (4.4) for C$ to uncover the universal entangle¬ 
ment entropy. Counter terms which are unavoidable even in the absence of interactions 
will contribute to (4.4) and result in a finite C^. This makes us believe that the univer¬ 
sal entanglement entropy for a QFT in the setup under study is directly related to the 
logarithmic divergences associated with renormalization of the bare operators and bare 
parameters in the energy-momentum tensor trace, whereas renormalization of the theory 
generates a finite EE through the use of (3.15). Of course, the finite part of EE depends 
on the choice of subtraction scheme, and therefore it is ambiguous, but the coefficients of 
the logarithmic running are universal. In what follows we elaborate the details of RG flow 
of the (finite) EE for a free massive scalar field in d = 2,4. These results generalize easily 
to higher dimensions. 

4.1.1 Renormalization of EE in even d 

Eq. (4.17) reveals how the universal EE is encoded in (3.15) (recall that dimensional 
regularization eliminates non-universal power law divergences). However, as we stressed 
in the previous section, this calculation is incomplete without accounting for the counter 
terms which are necessary to render the partition function of the theory finite. Taking 
these terms into account results in a finite vev of the energy-momentum tensor and thus 
finite EE through the use of (3.15). Of course, the universal divergence of EE is not lost. It 
transmutes into a finite term which represents the logarithmic running of EE. In this section 
we calculate all necessary counter terms needed to renormalize the partition function and 
obtain a renormalized expression for the EE. 

d = 2 

In d = 2, the total action which includes all possible counter terms that are necessary to 
render the partition function of the theory finite is 18 

StT 2 = J Q W) 2 + 7> m V + \^4> 2 + A 0 + K 0 llJ , (4.18) 

where Ao describes a bare cosmological constant and kq is the usual bare coupling of the 
Einstein action. In 2D this term is responsible for the trace anomaly. The conformal cou¬ 
pling £ c vanishes in two dimensions, however, we are going to use the method of dimensional 
renormalization, and therefore we keep it in the action. 

18 For a free field theory the mass and the non-minimal coupling are not corrected by the divergences of 
higher-order loops, therefore we do not distinguish between the bare and renormalized m 2 and £ c . 
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The bare and renormalized parameters in the minimal subtraction scheme are related 


as 

A » = ^ 2 ( A+ ^ 2 ” 2 ) ' 

k 0 = n d ~' 2 , ( 4 - 19 ) 

where A, k are renormalized parameters, // is an arbitrary mass scale and A p , k p are residues 
of simple poles. These counter terms are necessary to ensure a finite partition function. 
For an interacting QFT the above expression contains higher order poles and much more 
complicated residues. 

The constant A p is determined from the requirement that 

- m2 S = / s dT <02) + ^ 2m2 ^ E 7 <4 ' 20) 


must be finite. Expanding (4.9) around d = 2 gives 

= ~ 2tt (d- 2) + ~L ~ ^ “ A ) “ V>(A) - 7) + 0(d - 2) . (4.21) 


Hence, 


. 1 
A p — 

4-7T 


(4.22) 


Similarly, the free parameter k p is determined by the divergence of a 1-loop diagram 
with external gravitons. However, we resort to a different method. Let us consider the 
vacuum expectation value of the energy-momentum tensor trace 


(T(x)) = = _ ^Y “(^( _v2 + ^ + m 2 )^)-m 2 (0 2 )-dA o -(d-2)Ko^ . 

(4.23) 

The first term on the right hand side vanishes since it corresponds to the equation of motion 
operator’s vev. The rest can be expanded around d = 2 using (4.19), (4.21), (4.22) 


{T(x)) = - T ^— (log( 47 ri? 2 /x 2 ) - ^(1 - A) - 'ip(X) - 7 + l) - 2 A — k p TZ . (4.24) 

47T v 


As expected, the vev of the energy-momentum tensor trace is finite after the bare 
operators and bare parameters are replaced with renormalized parameters. The unspecified 
//, A and n p can be determined by imposing the decoupling condition ( T) —> 0 as mR 00 
[91] 


2 2 
[i = m 


A = m 


2 (7 - log( 47 r) - 1 ) 


87T 


K n = 


-1 

247T 


(4.25) 


Substituting into (T), we obtain 


<T(z)> 


9 

m 

47T 


(2 log(m-R) — -0(1 — A) — V’(^)) + 


n 

247T ■ 


(4.26) 


Note that the last term cannot be modified by adding to the action a finite counter term, 
hence it represents the trace anomaly in two dimensions. 
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Next we evaluate the EE flow using (3.8) and (3.15) 


C <t> = (2 log(mi?) - - 0(1 - A) - ip(X)) - ^ , 

dS EE ( m2 , *(™-R) 4 (^'(l-A)-^(A)) 1 

dR 2 yJ{mR) 2 - 1/4 3 


(4.27) 


It can be readily verified that RS' EE vanishes as mR —> oo. In this limit the RG running 
takes us to an empty fixed point. In the opposite limit, mR —> 0, we recover a pathological 
massless scalar field on S 2 (£ c = 0 in 2D) 19 . Therefore the above equations describe RG 
running of a massive theory into an empty fixed point. 


d = 4 

This case is similar to d = 2 except that the total action admits more counter terms, see 
e.g., [92], 

S?= 4 = J sd QW) 2 + \m 2 4> 2 + \tcKtf 

+ Ao + Ko 1Z + C^ upa Cp Upcr + 2a 0 E 4 + cqR 2 ^ , (4.28) 

where Cp Up(7 is the Weyl tensor. The relation between bare and renormalized parameters 
is 

A 0 = /-■> (a + Au nA , K 0 = M d - 4 (K+^m 2 ) , (4.29) 

a°=/- 4 (a+^), &0 = /-“(&+A_) . c° = K d -‘(c+jh) ■ 

As before, the constant A p is obtained by requiring that 

- m2 H=A (v ;,#, 2 )+ + 2 ^ 4 "* 4 ^i) (4 ' 30) 

is finite. Expanding (4.9) around d = 4 gives 

(0 2 > = 87 r 2 (d- 4 ) ~~ ^ ( 1 °g( 47r - R2 ) - V’(A) - V(3 - A) + 1 - 7 ) + 0(d - 4) . (4.31) 

Hence, 

Ap = ^ 9^2 ’ k p = 0 • (4.32) 

19 We call this theory ‘pathological’ since Laplace operator on a sphere is not invertible. It has a normal¬ 
izable zero mode which corresponds to a constant field configuration. Thus a minimally coupled massless 
scalar held on a sphere does not have a well-defined two point function without excluding zero mode from 
consideration. 
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To determine other counter terms, let us consider vev of the energy-momentum trace, 

(T(x)} = + i c 1l + m 2 )cp) - m 2 ((f) 2 ) - dA 0 - (d-2)n 0 TZ 

- (d - 4 ) (^2 C^ pa C^ pa + 2a 0 £ 4 + c 0 77 2 ^ + 4 (d - l)c 0 V 2 77 . (4.33) 

The first term is proportional to the equation of motion operator’s vev, thus it vanishes. 
Furthermore, 77 is constant on a sphere, whereas Weyl tensor vanishes. Thus, 


(T(x)) = —m 2 ((j) 2 ) - dA 0 - (d - 2 )k 0 K-(d- 4) (-2a 0 £ 4 + c 0 77 2 ) (4.34) 

4 Q 

= (log(4vri? 2 /r 2 ) - ip( A) - V’(3 - A) + - - 7 ) - 4A - 2kR - 2a p E 4 - c p H 2 , 


where we used (4.29), (4.31) and (4.32). The last two terms correspond to the trace 
anomaly. Obviously, these terms cannot be modified by adding a finite counter term to 
the action. 

Next we impose the decoupling condition, (T) —> 0 as mR —» 00 , to determine 
a p ,c p ,n,K, A [91] 


2 2 

H = m 


A = 


m 


4 3 - 2y + 21 og( 47 r) 
128vr 2 


K = 


—m 


(24tt) s 


27T 2 


Ob + 144 C„ = 


240t r 2 


(4.35) 


Although for our needs it is not necessary to calculate a p and c p separately, it is still 
worth mentioning that c p = 0 for a free massive scalar field. Indeed, c p is mass independent, 
and it vanishes for a massless conformally coupled scalar field (CFT). One can also verify 
that c p = 0 by a direct calculation [92, 93]. Substituting the above parameters into (T), 
we obtain 


(T(x)) = 


m 

167T 2 


2 log(mi?) — ip(X) — 3 — A) + 


3(mi ?) 5 


1 


240vr 2 i? 4 


(4.36) 


Finally, using (3.8) and (3.15), results in 


C, = - 


(mR) 4 

I 67 r 2 


21og(mi?) — ip(\) — V>(3 — A) + 


3 (mi ?) 5 


+ 


240vr 2 


R 


dS EE _ (mR) 2 (mR) 4 i(mR) 6 (tp'(3 - A) - ip'(\)) 


dR 


18 


12 


2AyJ(mR) 2 - 1/4 


1 

90 


(4.37) 


In the deep IR limit, mR — 00 , the theory flows into an empty fixed point and RS ' EE 
vanishes. In the UV limit, mR —> 0, we recover a conformally coupled massless scalar field, 
and RS ee = —1/90 in accordance with (3.32). 20 Therefore RG running happens between 
two fixed points which correspond to the conformally coupled scalar field and an empty 
theory. 

20 _ 

a — a p — 36Q . 
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4.2 Dirac fermion 

Let us consider a free Dirac field of mass m. The Euclidean action is given by 


I = d d x (y + m) 


and the corresponding energy-momentum tensor reads 

Tftv = - S a ^ytp + mipip) , 

where 7 ^ are the gamma matrices satisfying the anticommutation relations 

{ln,lv} = 2 g^u ■ 


(4.38) 


(4.39) 


(4.40) 


Taking the trace of the energy-momentum tensor and using the Dirac equation of motion 
we obtain, 

T = —rmfyip. (4-41) 


Hence 


(4.42) 


(4.43) 


= m R d (i/ji/j) . 

To evaluate C^ we use propagator of the Dirac field on S d [94] 

T(l + imR)F(l-imR) 

W( »» W> - 2 W(f + l)i^ 

X ^mRU(y, x) cos ^ 2 ^ 1 ^ +imR , ^ - imR ; ^ + 1 ; cos 2 

d .. TT , .. 0 ^ fd „ d ^ d o y\\ 

+ 2 n ‘ ^ U (y^ x ) sm 2 2jFi (,2 +imR ’ h _OTljR; 2 ’ cos 2 J j ’ 

where is the unit tangent vector to the geodesic connecting x to y, U(y,x ) is a matrix 
in the spinor indices which parallel propagates a spinor between the two points, and x is 
the polar angle between x and y. In particular, 


,7 M = [4] r (f + imR)T(l - imR) sinh(Trm) 

W i? rf_ 1 (47r ) rf / 2 sin(^ E )r (0 


(4.44) 


This vev is finite for odd d, whereas for even d it has simple poles which correspond to the 
logarithmic divergences. Let us consider these cases separately 

Odd d 

For odd d, we have 


d- 1 
2 


„ , . d+i (irmR) tanh( 7 rmii) -r-r ,. , . 2 n 2 \ 


>/2 (27r) d / 2 T ( 2 ) 


= H 


d+i (irmR) tanh( 7 rmf?) 
2 V2(27r) d / 2 T(0 


( ml ?)"- 1 + ^-^"^ ( mi ?)"- 3 


+ ...+ 


(4.45) 

r 2 (i)\ 


7T 
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The hyperbolic function is associated with the thermal corrections, and in principle it 
should be stripped off to isolate the impact of entanglement. We achieve this goal by 
taking the limit mR 3> 1, in which case the thermal effect is exponentially suppressed. To 
leading order we recover a known result [71], see also Appendix A 


a, = 


. d+i 

vr(-) 2 


A CDirac 
^CiOuniv 


«»i y/2{2n) d / 2 T (^) 


()T- 7r( d 2 ) 


(: mR) d + d{d ~ l ^ d ~ 2 \ mR) d ~ 2 


+ • 


r 2 (l) N 

. H- m R 


i r 


A T ,m d ~ 2 + ... , 


dR mR»i i2y/2(2n) — T (f) 
where ellipsis encapsulate curvature corrections. 


(4.46) 

(4.47) 


Even d 

The pole structure of for even do is given by 

rd 0 /2 


(27r)*/ 2 r (^ 


d — do 


+ • • 


(_)<fo/2-l 
(27r) rf o/2r (f 


( mR) a ° + 


dQ ^ do(do - l)(do - 2) dp-2 


24 


(4.48) 
2 


dp(d 0 - l)(d 0 - 2)(d 0 - 3)(d 0 - 4)(5d 0 + 2) rfo _ 4 + p 2 ( —) (mR) 2 \ _h 

5760 1 ’ ) d - do 

The universal entanglement entropy is obtained by substituting it into (3.15) 


R- 


clS'l 


irac 

univ 


(-) d o/2(d 0 - 2)A S / d 2 (d 0 — 3)(d 0 — 4)(5d 0 + 2) m d ° 4 


dR 6(27t) ( rf o—2)/2r 


m d °“" + 


120 


« 2 


+ ... + 12 


r 2 (f) 


m 


do (do - 1) R d ~ 4 


log 6 


(4.49) 


The leading order term matches a well-known universal ‘area law’ [71], see also Appendix A. 
Subleading terms represent corrections to the entanglement entropy induced by curvatures. 
In Appendix B we present independent computation of the curvature corrections and find 
full agreement with (4.49). Although we do not present it here, similar analysis to section 
4.1.1 for the massive scalar can be applied to the massive fermion, giving a renormalized 
expression for the EE, and the UV limit agrees with (3.32). 


5 Spectral decomposition 

In this section we derive an expression for the right hand side of (2.10) in terms of specific 
spectral function. We start from the spectral decomposition of the conserved two point 
function (3.10). In general, it can be expressed as [81] 

l(^/ii'(®)-^a/3(y))con ^'ofu/^afii'^i V) 4“ I/) > (5-1) 
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where the spin -2 piece, T 2 g,u,ai 3 , is traceless, and therefore it does not contribute to ( 2 . 10 ). 21 . 
The spin-0 piece is given by 

I'onu,ap(x, y) = S flu (x)F 0 (a) i S a ^(y) , 

S a /3 = Vq.V /3 - g a /3 V 2 - ^ - g a /3 , (5.2) 

where cr(x,y) is the geodesic interval in units of R between x and y , covariant derivatives 
on the left of ^(cr) act on x, whereas covariant derivatives on the right of To(er) act on y, 
as indicated by the arrow sign above S a p(y), and 

roc 

F 0 {a) = Q d -i / dgpo(y)Gfj,(a) , (5.3) 

where as before G m (ct) is the Green’s function satisfying 

( - V 2 + i c ll + y 2 )G^{x,y) = 5 d {x,y) , (5.4) 

with p being the mass of the held. Thus, the two point function which is necessary to 
evaluate the right hand side of ( 2 . 10 ) is completely determined by po(p-), 

fi;ti(T(x)T a/3 (y)) con = g^S^F^afs^y) . (5.5) 


Before substituting (5.5) into (2.10), it is worth mentioning that for our choice of 
foliation (2.4), we have £ /J ' = 5r and n^ = R cos 9 5^. Hence, 


V 2 = y TT V r V T + 


1 ^ / . d—2 x ^ 


r 2 gj n “ 2 q 


Z^dg (shl d 2 Odg) + 


R? sin 2 


-V 


S d-2 


i a n^S a p = R cos 0 ( 1 dg(sm d 2 Odg) + 

\ R 2 sin® 2 0 

with V ^ d _ 2 being the intrinsic Laplacian on S d ~ 2 . 

Next, we find it useful to define 


1 


(d- 1 ) 

R 2 sin 2 9 ' sd 2 " r R 2 


: V cd—2 + 


,(5.6) 


G [2) {t -t',0,6 1 ) = J dQ' G^t,0^t' ,0' ,&') (5.7) 

where (r, 9, fl) and (r 7 , O ', Vt') are coordinates (2.5) of x and y respectively, and the integral 
runs over a (d— 2)-dimensional sphere parametrized by f Y. The hrst thing to note about G ^ 
is that it depends on the difference r — t' due to the symmetry of (2.5) under translations 
in r, and it is independent of H since (2.5) is symmetric under rotations of S d ~ 2 . 
Combining (2.6) and (5.5), leads to 

(T(x)K) con = J + Jp + M 2 )po(/ 2 )G M (x,y)V a ^(y)^ a n /3 

= dfj,(£ c K+ +p 2 )po(p) (5.8) 

7T 

x J dO' cos 6' (dgi (sin d-2 O'dg /) + (d — 1 ) sin cZ ~ 2 0'^j G^ (r, 6, O') , 

21 Explicit expression for is derived in [81]. 
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where based on (5.4) we substituted V 2 G At (<r) = (£ c 7£ + g 2 )G fl (a) in the first equality, 22 
whereas in the second equality we used (5.6) and the definition (5.7). 23 We also assumed 
that d > 2. 24 Finally, integrating the right hand side of (5.8) by parts twice, yields 

(T(x)K) con = d^ c ll+ ^ + ^ 2 )po(^ 2) (T,M / = f) , (5-9) 

/ q\ 

where G^ { r, 9,9' = |) is independent of t since 9' = |, and the system exhibits rotational 
symmetry in the transverse space. As expected, (T(x)K) is just a function of 9 which 
parametrizes the geodesic distance from the entangling surface in the transverse space. 
Another useful result follows directly from (5.4) 


J d d xy/g{x)G^x,y) = - . (5.10) 

Indeed, the integral on the left hand side is convergent since the Green’s function is ev¬ 
erywhere regular on a sphere except at x = y, where it diverges as a 2 ~ d . However, this 
divergence is balanced by the integration measure which behaves as cA _1 . Furthermore, 
S d is a maximally symmetric space, and therefore the result of integration is given by some 
constant independent of y. In particular, (5.10) follows from the following identities 

(£ c K + y 2 ) J d d xy/g(x)G ll (x,y) = (-V 2 + £ C K + g 2 ) J d d x s/g{x)G^{x,y) = 1 , (5.11) 


where the last equality rests on (5.4). 

Substituting now (5.9) into (2.10), and using (5.7), (5.10) to integrate Glover x, 
yields 


R 


dS EE 

dR 


(d-i)r(g) 


d—2 
7T 2 


As / dg ( 1 + 




d 


R 2 ^ c n + g 2 ) 


Po(d) , 


(5.12) 


where As is the area of the entangling surface. 

The integral on the right hand side of (5.12) is convergent at the lower bound since 
sphere introduces a natural IR cut off, but it does not necessarily converges at the upper 
bound. Indeed, in the UV limit curvature corrections are irrelevant, and convergence of the 
integral is essentially the same as in Minkowski space. Note also that in the limit g?lZ ^$> 1 
we recover the result of [33]. 

Furthermore, if the integral on the right hand side converges, then positivity of the 
spectral function guarantees RS EE < 0 along the RG flow. However, it certainly does not 


22 We ignored the delta function on the right hand side of (5.4). There is nothing bad about it if x is 
disjoint from V. If, however, x hits the support of K, then one may question whether it is justified to 
suppress the delta function in (5.4). In general, if supports of the operators overlap, it is a must to include 
contact terms. However, inclusion of such contact terms will break the 0(2) invariance inherent to the 
setup under study, and therefore we disregard them. From computational point of view, we simply assume 
that x (j: V. 

23 Note that by assumption x is disjoint from V, therefore GR®, y)V| d _ 2 is everywhere regular on V and 
its integral over S d ~ 2 vanishes just because this manifold has no boundaries. 

24 The special case d = 2 can be treated similarly. In this case the entangling surface consists of two 
disjoint points 6' = ±^, and the integral over 9' in (5.8) should be extended from — ^ to ^ in order to 
cover all of S 2 , see (2.4). 
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ensure (4.15) since finite EE at the IR fixed point, 5^, includes terms which are part of 
the UV physics. These terms should be subtracted to isolate the universal contribution, 
Cq r . We elaborate details of this point in the discussion section. 


5.1 Example: conformally coupled scalar on S' 3 

In this section we carry out spectral decomposition for a massive conformally coupled scalar 
field described by the Euclidean action (4.1). In this case, based on (4.3), we have 

(T(x)T(y)) = m 4 (^ 2 (x)</> 2 (y)) = 2 m 4 G 2 n (x,y) , (5.13) 


where by assumption x and y are two disjoint points, and therefore the equation of motion 
operator in (4.3) does not contribute. In particular, it follows that the spin-0 piece of (5.1) 
is straightforwardly related to the spectral representation of 


( < l )2 (x)4> 2 (y)} = / dpp <j p(p)G fl (x,y) = dp? p 4> 2 (p)G fi {x,y) , (5.14) 

where p^i(p) = 2pp,p(p). In the meantime we keep the lower bound p^ unspecified. We 
find it convenient to use p 2 as the integration measure in the spectral representation since 
pp may admit imaginary values to account for the possibility of negative /i 2 on a sphere. 
A useful relation can be established between d m 2 (0 2 ) and p ^2 


- (4> 2 ) = J d d yyf^{(t> 2 {y)(t> 2 (x)) = j^ d P 2 ^^2 > 


dm 2 


(5.15) 


where in the last equality we used (5.10). Shifting the integration measure, 2 = ( Rp ) 2 — 
(Rpcp) 2 and using (4.9), yields 


'ip(X) — ip(d — 1 — A) — 7 r cot (f (d — 2 A)) 
M 


6 2 ) = iR 2 ~ d 


f 


dz 


p( z ) 


(d- l) 2 /4 + M 2 + z 


(5.16) 

where M 2 = (Rp^) 2 — \ and p ^2 = R' l ~ d p. We will use this result to check our computations 
in what follows. 

For the rest of this section we continue to explore the spectral decomposition in d = 3. 
In this case solution (4.7) to the Green’s equation can be expressed in terms of elementary 
functions 

sinh [M(x ~ 7r )] 


Gm(x) = - 


M=\l{mRY-\ 


(5.17) 


47 tR sinh( 7 rAA) sin y ’ 

We choose to work in the basis of scalar spherical harmonics on S 3 of radius R. They are 
given by [95] 


1 


R 3/2^ 

where Z 3 > I 2 > |ii| and 


lh4> 2 c\\P l2 h (cos 0) 3 cg (sin X ) 1/2 ^ 3 +i +a) ( cos x) , (5-18) 


PM(x) = 


iC L = 


1 


(£) 


M/2 


r(i + /i) 

2L + n- 1 (L + l + n — 2)H 1/2 

2 (L — l)\ 


2 F 1 \-v , Z 7 + 1 ; 1 + p ■ 


1 — x 


(5.19) 


- 24 - 
















Note, that by definition spherical harmonics satisfy 


Y Li^ e ^) = ^- 1 ) hY hMMx,0A) • 


(5.20) 


In particular, the only nontrivial coefficients in the expansion of G m (x), are given by 

R 1/2 h + 1 


Gh{m) = J s , Y,Mx)Gm(x) = V2*(h+ir +at 


where according to (5.18) 


Y h oo(x) = 


sin [(l 3 + l)x] 


i? 3 / 2 \/27r sin x 


Similarly, one can find the expansion of G^ix)- 
According to (5.17), 


G 2 m (x) = 


sinh 2 [A7(x — 7r)] 


167r 2 i? 2 sinh 2 (ttAT) sin 2 x 
Hence, its expansion in terms of spherical harmonics (5.18), is given by 

s j s3 YiMx)G 2 m {x) = ^7^^V^h^AO 

2AA sinh(27rAA) + • • • + 

Z§ + 4J\f 2 ) 

sinh {2tyJ\ f) f 1 


; l 3 odd 


— - 4 - 
2 ^ 


G\f(m 2 ) = 


2 7 / 2 R 1 / 2 tt 


1 — i 


. coth(-7rAA) r (/i , l 3 . K A , , h 


7r 


W 1 + tt - iN) - ip(l + -f + iN)\ 


(5.21) 


(5.22) 


(5.23) 


( 1 1 2A/* 1 2A/* 1 , 2A/ - '\ . ; 

4 + 1+A^ + 4+A^ h b (Ii72p+A^ J > ‘3 even 

Using the digamma function, t/>(a;), it can be succinctly written as 

1 


(5.24) 


(5.25) 


The spectral representation (5.14) is equivalent to the spectral representation in the 
angular momentum space 

1 /*°° 

G\f(m 2 ) = - di? p^{ji)Gi 3 {n 2 ) . (5.26) 

Substituting (5.21) and (5.25), and shifting the integration variable z = (Rfi) 2 — (i?//^) 2 , 
yields 


1 

4(Z 3 + 1) 


. coth(7rAT) 

1 — 2- 

7r 


Wi + | 


W) - -0(1 + + W)] 



p(z) 

(h + l) + A/"| + z 


(5.27) 


As expected, the special case of this formula, l 3 = 0, matches (5.16). 

Before we proceed, let us discuss possible values of /i 2 . If /r 2 = 0 it means that the 
spectrum of particles starts from a conformally coupled scalar field (z = (1 in this case) and 
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continues all the way up to infinite massive modes (z = oo). However, there is nothing bad 
about negative /r 2 since conformal coupling in the Euclidean action (4.1) may compensate 
its negativity such that the overall q i 2 term is positive. The lowest possible negative value 
of /r 2 is given by (Rp ^) 2 = — It can be read off either from the Euclidean action (4.1), 
or by setting z = I 3 = 0 in the integrand of (5.27) and demanding positivity. 

We need to invert (5.27) to get the spectral density p(z). For brevity, we define 


zo — (h + l) + -A/^ =>■ l 3 --l + y/zo- A/"|, 

/(«.)= 

Jo z 0 + Z 

To express p(z) in terms of f(z), we first note that by definition 
v f(-zo-ie)-f(-zo + ie) 


lim 

e-s-0 

Combining this result with 


2 i 


= lim / dz p(z) 


£—^0 


(z - zo ) 2 + e 2 


S(z — Zo) = - lim 7 - %- , 

v ' 7r e->o (z - z 0 ) 2 + e 2 


(5.28) 


(5.29) 


(5.30) 


yields, 


iHz) = 1 lim fi-z-V-n-t + K) 
' 7T £— >0 2i 


Using now (5.27), (5.28) and the definition of p, we get, 


(5.31) 


pR f ^ . coth(7rA/") 

~ 2 lT 




^(^(! - iMu) - i/V) - V>(^(1 - itfn) + iN) 


+V>(^(1 + iJJ'n) ~ iN) ~ ^(^(1 + iNfi) + iN) 


®((pR) 2 -l/A) (5.32) 


where Af tJ = -y/ (pR ) 2 — 1/4. There are certain limits when (5.32) can be checked. 


In the limit of conformally coupled scalar which is also a UV limit, mR —> 0, and we 
have J\T = |. Hence, in this limit 


Pqp(fJ) 


pR 


-®{JpR ) 2 — 1/4) . 


mR-^o 2-K^{pR ) 2 - 1/4 
and it can be checked by a direct computation that (5.27) holds. 

Flat space limit is recovered if we fix m and p while R —> 00 , i.e., R 3> m -1 , p~ l 


(5.33) 


o 27T 


(5.34) 


where we used ip(x) ~ log(x) for i»l. In particular, in accord with [96], we obtain 


Mf,) rIoo 


(5.35) 
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Substituting into (5.12) gives RS EE ~ mR S> 1. On the other hand, according to 
(2.10) one expects RS' EE = 0 since m.R —> oo limit represents deep IR in our setup, 
which corresponds to an empty theory in this case. This contradiction is, however, 
apparent. Indeed, R represents the characteristic RG scale and therefore one should 
only consider /r ~ 1/R <C m. In this region theta-function and hence po vanish. On 
the other hand, the spectral integral (5.12) accounts for all possible scales including 
/r m, and therefore RS EE is contaminated by various terms which are not part of 
the IR physics. 


• It worth noting that in the limit of fiR 3> 1 with m and R fixed, we expect that the 
spectral function on a sphere asymptotically approaches its counterpart in flat space. 
Indeed, 


P<j> *(/*) 


1 

1 27T 


(5.36) 


6 Discussion 


In this paper we study the renormalization group flow of entanglement entropy in an an¬ 
alytically tractable setup - a cap-like entangling region in de Sitter space. Our system 
has much in common with the indispensable Rindler wedge, e.g., the entangling surface 
exhibits 0(2) symmetry in the transverse space, and the corresponding entanglement en¬ 
tropy equals thermal entropy for any QFT. However, in contrast to the Rindler space, 
where the geometry is flat and there is no characteristic temperature, in our setup the 
background curvature does not vanish, and it sets a characteristic length scale, R , for both 
the curvature of the entangling surface and temperature of the environment. 

Since there is only one global geometric scale, it also effectively determines the mass 
scale of RG running to be of order 1/R. In particular, studies of RG flow boil down to a 
constant Weyl rescaling of the background geometry. From this perspective, our setup is 
scalable [26], and we argue that EE satisfies the following RG equation 


dSvE _ V S d d 
R ^m ~ T R dR^' 


( 6 . 1 ) 


with Vgd being the volume of a d-dimensional sphere. 

The above simple relation between the entanglement entropy flow and trace of the 
energy-momentum tensor allows to analytically explore various properties of the entangle¬ 
ment entropy running when the system flows between the UV and IR fixed points. 

In section 4 we scrutinize the RG flows of EE in d = 2,3,4 for a conformally coupled 
scalar field deformed by a mass operator. In d = 3 our findings indicate that the renormal¬ 
ized entanglement entropy defined in [26, 28, 54] exhibits stationarity at the fixed points 
on a sphere, but it is not monotonic along the RG trajectory. This result is completely 
opposite to the behaviour of REE in Minkowski space, where it has been shown that REE 
is monotonic [54], but not stationary [69]. 
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Lack of monotonicity may result from various causes. For instance, REE in de Sitter 
space receives thermal corrections which are absent in Minkowski space. Furthermore, ab¬ 
sence of IR divergences in a QFT living on a sphere might be another plausible explanation 
for a qualitative difference in the behaviour of REE on de Sitter and on Minkowski space. 
That being said these arguments are unlikely to explain the discrepancies in the behaviour 
of REE in the vicimity of UV fixed point, since in the deep UV both the curvature and 
temperature of de Sitter space are irrelevant. 

Moreover, our calculations reveal another interesting and closely related result. We 
find that the universal ‘area law’ for a conformally coupled scalar field is different from the 
known expression in Minkowski space [71-73]. This discrepancy cannot be attributed to 
a curved geometry since the universal ‘area law’ depends solely on the area of entangling 
surface. This observation raises a controversial [15, 33, 36, 74-78] question whether there 
exists any difference between the universal entanglement entropies for minimally and con¬ 
formally coupled scalar fields in flat space. Possible explanation and interpretation of this 
phenomenon can be found in [79] and [56, 62, 80], see also [33] for interpretation based on 
calculations of the mutual information in Minkowski space. Our findings here point out in 
favour of the difference between the minimally and conformally coupled scalar fields in the 
limit of flat space. 

In section 5 we derive the spectral representation of entanglement entropy flow. The 
final expression (5.12) depends on a specific spectral function which determines the spin-0 
part of a two-point function of the energy-momentum tensor [81]. Reflection positivity 
ensures that the spectral function is positive, and therefore (5.12) suggests that the rate 
of change of entanglement entropy along the RG trajectory is negative provided that the 
integral on the right hand side converges. In particular, it implies > S'™ for a finite 
part of EE in d = 3. 

Obviously, this inequality is not the same as (4.15), and therefore it does not prove 
the E-theorem [52, 53] in three dimensions. To understand the difference between the two 
inequalities, it is enough to consider the universal ‘area law’ for a massive free field, see 
Appendix A. In d = 3 it represents a finite contribution to EE which grows linearly in the 
IR limit. Since massive degrees of freedom decouple in the deep IR, the universal ‘area law’ 
becomes part of the UV physics, and therefore it should be subtracted to isolate a ‘true’ 
universal entanglement entropy in this limit. 

Of course, IR entanglement entropy may contain all kind of finite terms associated with 
the fingerprints of UV physics. Such remnants of the RG trajectory should be subtracted 
to extract the universal piece. REE is a particular subtraction scheme which proved to be 
effective in building a c-function in Minkowski space [54]. However, we have shown that 
it is not as good on a sphere. Unfortunately, we were not able to identify a reasonable 
candidate for a c-function on a sphere. Thermodynamic inequalities might be a good source 
to look for plausible candidates, e.g., the generalized second law is one such example [97]. 
We will explore this avenue elsewhere. 
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A Universal area law 

In has been shown in [33, 62] that variation of the entanglement entropy with respect to 
the relevant coupling can be expressed in terms of a particular structure in the spectral 
decomposition of the energy-momentum tensor 



(A.l) 


where (/r) is the spectral function which corresponds to spin s = 0 states in Minkowski 
space [96]. For free conformally coupled scalar and Dirac fermion, we have [96] 



0(/2 — 2m) . 


@(/r — 2m) , 


(A.2) 


Substituting into (A.l), and using dimensional regularization, yields 



(A.3) 



(A.4) 


Hence, we get 




d 


U- — 4 / J \ 

6 ( 2 tt) T-r(§) 

m = < 



is m d 2 log(m<5) for even d , 



(A.5) 


k 12V2(27T)^r(f) 

where simple pole in the gamma function is identified with log(md) 25 
Similarly for the conformally coupled scalar, 



for even d 


for odd d. 


(A.6) 


25 To establish this dictionary, it is enough to introduce a sharp cut off p, 


'max 


1/5 in (A.l). 
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Note that this result is different from the universal ‘area law’ for minimally coupled scalar 
field [71-73], 


rn 2 dS »^ = 

dm 2 


-As m d 2 log(m<5) 

12 (47r)Vr(f) V ' 

jd-2) , d -2 

24 d—2 . rn 

(4?r)U3-r(f) 


for even d , 
for odd d. 


(A.7) 


B Free fields on a deformed waveguide 

In this Appendix we perform an independent computation of the leading curvature cor¬ 
rection displayed in (4.17) and (4.49). To this end one can resort to the replica trick 
[7, 11, 98, 99] combined with the heat kernel technique, see e.g., [100]. Indeed, rotational 
symmetry around the entangling surface makes it possible to introduce a well-defined con¬ 
ical defect [51, 70, 78]. However, we choose to follow a different approach [61], where a 
special role played by the neighbourhood of the entangling surface is made manifest. 

Recall that divergences of the entanglement entropy are local and dominated by the 
region near the entangling surface. In particular, to recover (4.17) and (4.49) it is enough 
to explore our metric (2.5) in the vicinity of the entangling surface at 6 = it /2 

2 2 

ds 2 = (1 - ^2 + • • O 7 " 2 ^ r2 + R 2 d0 2 + (1 - -^2 + ... )R 2 d£l 2 d _ 2 , (B.l) 

where r = R{9 — vr/2). To leading order in the radial distance r, the geometry can be 
approximated by a waveguide with spherical cross-section R 2 x S d ~ 2 , whereas subleading 
terms can be treated as perturbations. 

In Cartesian coordinates x 1 = rcos9 ,x 2 = rsinO, the expansion of the metric takes 
the form 

ds 2 = (1 - ^RabcdX b x d + .. .)dx a dx c + ( 7 + n iacj x a x c + .. .)dy l dy j , (B.2) 

where ellipsis encode higher derivative terms which are irrelevant for our needs, {y l } d Z 2 
and 7 ij are coordinates and induced metric on S d ~ 2 respectively, and IZ^uap is the Riemann 
curvature tensor on S d ^ 

'R'p.vpo = ~^2 {dp.pdi'tr 9pcr9up) ■ (B.3) 

Of course, for a generic entangling surface (B.2) contains terms with extrinsic curvatures 
[16, 61, 101-104], 26 however in our setup they vanish due to rotational symmetry around 
S. 

Next we observe that the curvature terms in (B.2) are small close to the entangling 
surface (x a = 0 ) where the divergences are localized, and therefore we treat them as 
perturbations, h^, of the waveguide geometry, see [61] for details. In particular, it follows 
from ( 2 . 8 ) that linear correction to the entanglement entropy takes the form 

6 S ee = \Jd 2 x J d d - 2 y^(T^(x,y)K)h^(x,y) + 0 {h 2 ) . (B.4) 

26 See also critique of [103, 104] in [105]. 


- 30 








Only universal divergences of the above expression can be taken at face value, otherwise 
there is no reason to expect that higher order terms in h^ u do not contribute. 

Now the first term within parenthesis in (4.17), (4.49) represents the universal ‘area 
law’ [71-73], see also [33, 62] and Appendix A. This term is entirely encoded in the lead¬ 
ing order entanglement entropy, whereas SS EE contributes to the second term in (4.17), 
(4.49) which is proportional to the Riemann tensor. Given that h^ u is also proportional to 
the Riemann tensor, we can replace (T ,llJ K) by its flat space counterpart without losing 
contributions to (4.17), (4.49), then according to [62] 


6 S™ = 


7 T 12^—1 


2 d ~ 2 (d - l) 2 {d + l)T(d) 


S ac 6 ^n iajc + - 6 ° 


s bd n abcd 


% W (d) (B.5) 
d 2 


+ 




2 d ~ 2 {d- l)(d+ l)T(d) 


d-2 


s ac s bd n abcd - s ac s i] n. 


icljc 


dfi o(2) 


/ o d 




where and c^ 2 \^l) are two spectral functions which define a two point function of 

the energy-momentum tensor for a generic QFT in Minkowski space [96]. For free fields 
c(°)(//) is given by (A.2), whereas c^ 2 \^l) reads [96] 


c <?( k ) = 2 | i / 21 7 ^/- 3 (l 


4 m 2 
d 2 


2 \ (d~ l )/2 


2 4m 2 ', 

1 + ~ 2m) • 


<fV) = 


Hence, 

5 S' 


d 


d -3 


1 


4r7)2 x (d+m 

1 - ) 0 (/i — 2 m) . 


d 


(B.6) 


'scalar 

univ 


(-)V(d- 2)(-76 + 84d-25d 2 + 2d 3 ) A E m' 


d—4 


240(47r)^((2- l) 2 r(^) 


R 2 


log(md) , 


6 S, 


Dirac 

univ 


(-)5(d-2)(d-3) A s m d ~ 4 

60 ( 2 vr)^r(f) 


log(m<5) , 


(B.7) 


where we introduced a UV cut off /i max ~ 1/5. This contribution should be combined 
with the leading order entanglement entropy. To leading order the geometry is identical to 
a waveguide with spherical cross-section, and therefore for Dirac fermion we can use the 
results of [73] 




Dirac 

univ 


✓ x a 

(-1)2 A s 


R 2 xS d ~ 2 


6(2vr) 


d-2 

2 


ni 


m d -2 + (d ~ 2 )( d ~ 3 ) m> 


d—4 


24 


R 2 


+ • • • ) log(m<5) , (B.8) 


while for conformally coupled scalar, the result reads 2 ' 


S' 


scalar 

univ 


R 2 xS d ~ 2 


(- 1 ) 2 ( 1-6 g c )A s 

6(47r)^r(f) 


m d ~ 2 — 


(d — 2) 2 {d — 3)(1 — 6£ c ) m 


d—4 


12 


R 2 


+ .. 


log(m<5) , 
(B.ll) 


27 Entanglement entropy of a scalar field on a waveguide geometry was studied in [73] using the heat 
kernel methods. In particular, the authors apply the same heat kernel, Kc n , on a two-dimensional cone, C„, 
with an angular excess 2n(n — 1) to evaluate entanglement entropy for both minimally and non-minimally 
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Combining altogether, we obtain 


5 


'scalar 

univ 


5 


'scalar 

univ 


R 2 xS d ~ 2 


+ <55 


scalar 

univ 


(do - 4) (-^o/ 2 +l^ 

12 (d - 1 ) ( 4 7 r )(d 0 - 2 )/ 2 p f<k 


x 



(do - 2)(d 0 - 6 ) (5 d? 0 - 18d 0 + 4) m d °- 4 \ 
120(d — 4) W~ J 


log (md) , 


S Dirac QDirac 


R 2 xS d ~ 2 


+ <55, 


Dirac 

univ 


(B.12) 


(-1)U E / d _ 2 , (d-2)(d-3)(5d + 2) m d - 4 n , „ 

^ + 11 —^r- + ... log(m<5) , 


6(2vr)^r(|) 


120 


R 2 


in full agreement with (4.17), (4.49). 

C Addition theorem for spherical harmonics on S 3 


In this appendix we review the proof of the following useful mathematical result 

>boo( 7 ) = rfiE E KhiM Y hhh(x) , (c.i) 


3+ l 2 =0h=-l 2 

where x,y are two arbitrary points on 5 3 and 7 is an angle between them. This identity 
is called the addition theorem for spherical harmonics. 

Any function f(x) on 5 3 can be expanded as follows 


OO Z3 I2 

f ( x ) = E E E 

Z3—0 Z2 —0 h=—l2 


= / Khu ( x )f( x ) 


(C.2) 


If we set x = 0 then according to (5.18),(5.19) only terms with l 2 = 0 => l\ = 0 contribute 

fa + 1 


/(x)lx=o = E ^300 

h=o 


y/ 2 i t 


Ai 3 00 = / Yl 3 oo(x)f{x) ■ 


(C.3) 


On the other hand, expanding Y) 3 oo ( 7 ) in spherical harmonics and taking into account 
that it is a spherical harmonic of order fa, yields 

h h f 

^300(7) = E E ■^hhh (y)Yl 3 l 2 h O^) > Ahhh id) I hoo(l') ■ (C- 4 ) 

, O , , j X 


h =0 h=—fo 


(B.9) 


coupled scalar fields. However, to get (B.ll) we make use of the following expression 
TrK Cn (t) = (1 ~ n) (l - 60 + nTrA +2 (t) + 0(1 - nf , 

where £ represents a non-minimal coupling to the background curvature. To derive it, there is no need to 
find or use the exact solution [106-108] to the heat kernel equation on a cone. Instead, one can use the 
general form of the heat kernel expansion in the vicinity of t = 0 [100], and substitute curvature scalar on 
C n 

7Z = 4tt(1 - n)<5 s + 0(1 - n ) 2 , (B.10) 

where fe is a two-dimensional 5-function supported on E. 
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A hhh {y) can be viewed as Ai 3 oo coefficient in an expansion of Y*i 2 li {x) in a series of 
^ 3 Mi ( 7 ) referred to the axis y. From (C.3) it is then found that, since only one I 3 is 
present this coefficient is given by 

A hhiAy) = j^i Y hi 2 h(y) > ( c - 5 ) 

from which (C.l) follows. 

References 

[1] A. Kitaev and J. Preskill, “Topological entanglement entropy,” Phys.Rev.Lett. 96 (2006) 
110404, arXiv:hep-th/0510092 [hep-th]. 

[2] M. Levin and X.-G. Wen, “Detecting topological order in a ground state wave function,” 
Phys. Rev. Lett. 96 (Mar, 2006) 110405. 

http://link.aps.org/doi/10.1103/PhysRevLett.96.110405. 

[3] B. Hsu, M. Mulligan, E. Fradkin, and E.-A. Kim, “Universal entanglement entropy in 
two-dimensional conformal quantum critical points,” Phys. Rev. B 79 (Mar, 2009) 115421. 
http://link.aps.org/doi/10.1103/PhysRevB.79.115421. 

[4] T. Grover, A. M. Turner, and A. Vishwanath, “Entanglement Entropy of Gapped Phases 
and Topological Order in Three dimensions,” Phys.Rev. B84 (2011) 195120, 

arXiv:1108.4038 [cond-mat.str-el] . 

[5] M. Srednicki, “Entropy and area,” Phys.Rev.Lett. 71 (1993) 666-669, 
arXiv:hep-th/9303048 [hep-th]. 

[6] C. Holzhey, F. Larsen, and F. Wilczek, “Geometric and renormalized entropy in conformal 
field theory,” Nucl.Phys. B424 (1994) 443-467, arXiv:hep-th/9403108 [hep-th]. 

[7] P. Calabrese and J. L. Cardy, “Entanglement entropy and quantum field theory,” 
J.Stat.Mech. 0406 (2004) P06002, arXiv:hep-th/0405152 [hep-th]. 

[8] H. Casini and M. Huerta, “Entanglement entropy in free quantum field theory,” J.Phys. 
A42 (2009) 504007, arXiv:0905.2562 [hep-th], 

[9] P. Calabrese and J. Cardy, “Entanglement entropy and conformal field theory,” J.Phys. 
A42 (2009) 504005, arXiv:0905.4013 [cond-mat.stat-mech] . 

[10] L. Bombelli, R. K. Koul, J. Lee, and R. D. Sorkin, “A Quantum Source of Entropy for 
Black Holes,” Phys.Rev. D34 (1986) 373-383. 

[11] J. Callan, Curtis G. and F. Wilczek, “On geometric entropy,” Phys.Lett. B333 (1994) 
55-61, arXiv:hep-th/9401072 [hep-th]. 

[12] S. Ryu and T. Takayanagi, “Holographic derivation of entanglement entropy from 
AdS/CFT,” Phys.Rev.Lett. 96 (2006) 181602, arXiv:hep-th/0603001 [hep-th]. 

[13] S. Ryu and T. Takayanagi, “Aspects of Holographic Entanglement Entropy,” JHEP 0608 
(2006) 045, arXiv:hep-th/0605073 [hep-th]. 

[14] V. E. Hubeny, M. Rangamani, and T. Takayanagi, “A Covariant holographic entanglement 
entropy proposal,” JHEP 0707 (2007) 062, arXiv:0705.0016 [hep-th], 

[15] S. N. Solodukhin, “Entanglement entropy of black holes,” Living Rev.Rel. 14 (2011) 8, 
arXiv:1104.3712 [hep-th]. 


- 33 - 



[16] A. Lewkowycz and J. Maldacena, “Generalized gravitational entropy,” JHEP 1308 (2013) 
090, arXiv:1304.4926 [hep-th] . 

[17] T. Faulkner, A. Lewkowycz, and J. Maldacena, “Quantum corrections to holographic 
entanglement entropy,” JHEP 1311 (2013) 074, arXiv: 1307.2892. 

[18] T. Faulkner, M. Guica, T. Hartman, R. C. Myers, and M. Van Raamsdonk, “Gravitation 
from Entanglement in Holographic CFTs,” JHEP 1403 (2014) 051, arXiv: 1312.7856 
[hep-th] . 

[19] R. Bousso, H. Casini, Z. Fisher, and J. Maldacena, “Proof of a Quantum Bousso Bound,” 
Phys.Rev. D90 no. 4, (2014) 044002, arXiv: 1404.5635 [hep-th]. 

[20] O. Ben-Ami, D. Carmi, and J. Sonnenschein, “Holographic Entanglement Entropy of 
Multiple Strips,” JHEP 1411 (2014) 144, arXiv: 1409.6305 [hep-th]. 

[21] N. Engelhardt and A. C. Wall, “Quantum Extremal Surfaces: Holographic Entanglement 
Entropy beyond the Classical Regime,” JHEP 1501 (2015) 073, arXiv: 1408.3203 
[hep-th] . 

[22] H. Casini and M. Huerta, “A Finite entanglement entropy and the c-theorem,” Phys.Lett. 
B600 (2004)142-150, arXiv:hep-th/0405111 [hep-th]. 

[23] R. C. Myers and A. Sinha, “Seeing a c-theorem with holography,” Phys.Rev. D82 (2010) 
046006, arXiv:1006.1263 [hep-th]. 

[24] R. C. Myers and A. Sinha, “Holographic c-theorems in arbitrary dimensions,” JHEP 1101 
(2011)125, arXiv:1011.5819 [hep-th], 

[25] R. C. Myers and A. Singh, “Comments on Holographic Entanglement Entropy and RG 
Flows,” JHEP 1204 (2012) 122, arXiv: 1202.2068 [hep-th], 

[26] H. Liu and M. Mezei, “A Refinement of entanglement entropy and the number of degrees of 
freedom,” JHEP 1304 (2013) 162, arXiv: 1202.2070 [hep-th]. 

[27] S. N. Solodukhin, “The a-theorem and entanglement entropy,” arXiv: 1304.4411 
[hep-th] . 

[28] H. Liu and M. Mezei, “Probing renormalization group flows using entanglement entropy,” 
JHEP 1401 (2014) 098, arXiv: 1309.6935 [hep-th], 

[29] S. Banerjee, A. Bhattacharyya, A. Kaviraj, K. Sen, and A. Sinha, “Constraining gravity 
using entanglement in AdS/CFT,” JHEP 1405 (2014) 029, arXiv: 1401.5089 [hep-th], 

[30] S. N. Solodukhin, “Conformal a-charge, correlation functions and conical defects,” 

Phys.Lett. B736 (2014) 283-287, arXiv: 1406.5368 [hep-th]. 

[31] T. Kawano, Y. Nakaguchi, and T. Nishioka, “Holographic Interpolation between a and F” 
JHEP 1412 (2014) 161, arXiv: 1410.5973 [hep-th], 

[32] V. Balasubramanian, J. J. Heckman, and A. Maloney, “Relative Entropy and Proximity of 
Quantum Field Theories,” arXiv: 1410.6809 [hep-th], 

[33] H. Casini, F. Mazzitelli, and E. T. Lino, “Area terms in entanglement entropy,” 
arXiv:1412.6522 [hep-th], 

[34] M. Nozaki, T. Numasawa, and T. Takayanagi, “Quantum Entanglement of Local Operators 
in Conformal Field Theories,” Phys.Rev.Lett. 112 (2014) 111602, arXiv: 1401.0539 
[hep-th]. 


- 34 - 


[35] J. Cardy and C. P. Herzog, “Universal Thermal Corrections to Single Interval 
Entanglement Entropy for Two Dimensional Conformal Field Theories,” Phys.Rev.Lett. 

112 no. 17, (2014) 171603, arXiv:1403.0578 [hep-th] . 

[36] C. P. Herzog, “Universal Thermal Corrections to Entanglement Entropy for Conformal 
Field Theories on Spheres,” JHEP 1410 (2014) 28, arXiv: 1407.1358 [hep-th] . 

[37] L.-Y. Hung, R. C. Myers, and M. Smolkin, “Twist operators in higher dimensions,” JHEP 
1410 (2014) 178, arXiv:1407.6429 [hep-th], 

[38] T. Faulkner, “Bulk Emergence and the RG Flow of Entanglement Entropy,” 
arXiv:1412.5648 [hep-th]. 

[39] M. Goykhman, “On entanglement entropy in’t Hooft model,” arXiv: 1501.07590 
[hep-th] . 

[40] X. Huang, L.-Y. Hung, and F.-L. Lin, “OPE of the stress tensors and surface operators,” 
arXiv:1502.02487 [hep-th]. 

[41] R.-X. Miao, “Universal Terms of Entanglement Entropy for 6d CFTs,” arXiv: 1503.05538 
[hep-th] . 

[42] A. Zamolodchikov, “Irreversibility of the Flux of the Renormalization Group in a 2D Field 
Theory,” JETP Lett. 43 (1986) 730-732. 

[43] J. L. Cardy, “Is There a c Theorem in Four-Dimensions?,” Phys.Lett. B215 (1988) 749-752. 

[44] Z. Komargodski and A. Schwimmer, “On Renormalization Group Flows in Four 
Dimensions,” JHEP 1112 (2011) 099, arXiv: 1107.3987 [hep-th]. 

[45] Z. Komargodski, “The Constraints of Conformal Symmetry on RG Flows,” JHEP 1207 
(2012) 069, arXiv:1112.4538 [hep-th], 

[46] M. A. Luty, J. Polchinski, and R. Rattazzi, “The a-theorem and the Asymptotics of 4D 
Quantum Field Theory,” JHEP 1301 (2013) 152, arXiv: 1204.5221 [hep-th]. 

[47] H. Elvang, D. Z. Freedman, L.-Y. Hung, M. Kiermaier, R. C. Myers, et al., “On 
renormalization group flows and the a-theorem in 6d,” JHEP 1210 (2012) Oil, 
arXiv:1205.3994 [hep-th]. 

[48] S. N. Solodukhin, “Entanglement entropy, conformal invariance and extrinsic geometry,” 
Phys.Lett. B665 (2008) 305-309, arXiv:0802.3117 [hep-th]. 

[49] H. Casini and M. Huerta, “Entanglement entropy for a Maxwell field: Numerical 
calculation on a two dimensional lattice,” arXiv: 1406.2991 [hep-th] . 

[50] H. Casini, M. Huerta, and R. C. Myers, “Towards a derivation of holographic entanglement 
entropy,” JHEP 1105 (2011) 036, arXiv: 1102.0440 [hep-th]. 

[51] J. Dowker, “Entanglement entropy for odd spheres,” arXiv: 1012.1548 [hep-th]. 

[52] D. L. Jafferis, I. R. Klebanov, S. S. Pufu, and B. R. Safdi, “Towards the F-Theorem: N=2 
Field Theories on the Three-Sphere,” JHEP 1106 (2011) 102, arXiv: 1103.1181 [hep-th], 

[53] I. R. Klebanov, S. S. Pufu, and B. R. Safdi, “F-Theorem without Supersymmetry,” JHEP 
1110 (2011) 038, arXiv:1105.4598 [hep-th], 

[54] H. Casini and M. Huerta, “On the RG running of the entanglement entropy of a circle,” 
Phys.Rev. D85 (2012) 125016, arXiv: 1202.5650 [hep-th]. 


- 35 - 


[55] E. Lieb and M. Ruskai, “Proof of the strong subadditivity of quantum-mechanical entropy,” 
J.Math.Phys. 14 (1973) 1938-1941. 

[56] V. Rosenhaus and M. Smolkin, “Entanglement Entropy for Relevant and Geometric 
Perturbations,” JHEP 1502 (2015) 015, arXiv: 1410.6530 [hep-th] . 

[57] V. Rosenhaus and M. Smolkin, “Entanglement Entropy Flow and the Ward Identity,” 

Phys.Rev.Lett. 113 no. 26, (2014) 261602, arXiv: 1406.2716 [hep-th]. 

[58] J. J. Bisognano and E. H. Wichmann, “On the duality condition for a Hermitian scalar 
field,” Journal of Mathematical Physics 16 no. 4, (1975) . 

[59] J. J. Bisognano and E. H. Wichmann, “On the duality condition for quantum fields,” 
Journal of Mathematical Physics 17 no. 3, (1976) . 

[60] D. N. Rabat and M. Strassler, “A Comment on entropy and area,” Phys.Lett. B329 (1994) 
46-52, arXiv:hep-th/9401125 [hep-th]. 

[61] V. Rosenhaus and M. Smolkin, “Entanglement Entropy: A Perturbative Calculation,” 
JHEP 1412 (2014) 179, arXiv: 1403.3733 [hep-th], 

[62] V. Rosenhaus and M. Smolkin, “Entanglement entropy, planar surfaces, and spectral 
functions,” JHEP 1409 (2014) 119, arXiv: 1407.2891 [hep-th], 

[63] S. Banerjee, “Wess-Zumino Consistency Condition for Entanglement Entropy,” 

Phys.Rev.Lett. 109 (2012) 010402, arXiv: 1109.5672 [hep-th]. 

[64] A. Allais and M. Mezei, “Some results on the shape dependence of entanglement and Renyi 
entropies,” arXiv: 1407.7249 [hep-th] . 

[65] M. Mezei, “Entanglement entropy across a deformed sphere,” arXiv: 1411.7011 [hep-th]. 

[66] A. Lewkowycz and E. Perlmutter, “Universality in the geometric dependence of Renyi 
entropy,” arXiv:1407.8171 [hep-th]. 

[67] J. Maldacena and G. L. Pimentel, “Entanglement entropy in de Sitter space,” JHEP 1302 
(2013) 038, arXiv:1210.7244 [hep-th], 

[68] W. Fischler, S. Kundu, and J. F. Pedraza, “Entanglement and out-of-equilibrium dynamics 
in holographic models of de Sitter QFTs,” JHEP 1407 (2014) 021, arXiv: 1311.5519 
[hep-th]. 

[69] I. R. Klebanov, T. Nishioka, S. S. Pufu, and B. R. Safdi, “Is Renormalized Entanglement 
Entropy Stationary at RG Fixed Points?,” JHEP 1210 (2012) 058, arXiv: 1207.3360 
[hep-th] . 

[70] J. Dowker, “Sphere Renyi entropies,” J.Phys. A46 (2013) 225401, arXiv: 1212.2098 
[hep-th] . 

[71] M. P. Hertzberg and F. Wilczek, “Some Calculable Contributions to Entanglement 
Entropy,” Phys.Rev.Lett. 106 (2011) 050404, arXiv: 1007.0993 [hep-th]. 

[72] M. Huerta, “Numerical Determination of the Entanglement Entropy for Free Fields in the 
Cylinder,” Phys.Lett. B710 (2012) 691-696, arXiv: 1112.1277 [hep-th]. 

[73] A. Lewkowycz, R. C. Myers, and M. Smolkin, “Observations on entanglement entropy in 
massive QFT’s,” JHEP 1304 (2013) 017, arXiv: 1210.6858 [hep-th]. 

[74] S. N. Solodukhin, “Nonminimal coupling and quantum entropy of black hole,” Phys.Rev. 
D56 (1997) 4968-4974, arXiv:hep-th/9612061 [hep-th]. 


- 36 - 


[75] M. Hotta, T. Kato, and K. Nagata, “A Comment on geometric entropy and conical space,” 
Class. Quant. Grav. 14 (1997) 1917-1925, arXiv:gr-qc/9611058 [gr-qc] . 

[76] T. Nishioka, “Relevant Perturbation of Entanglement Entropy and Stationarity,” Phys.Rev. 
D90 (2014) 045006, arXiv: 1405.3650 [hep-th]. 

[77] J. Lee, A. Lewkowycz, E. Perlmutter, and B. R. Safdi, “Renyi entropy, stationarity, and 
entanglement of the conformal scalar,” arXiv: 1407.7816 [hep-th] . 

[78] J. Dowker, “Expansion of Renyi entropy for free scalar fields,” arXiv: 1408.4055 [hep-th]. 

[79] R. C. Myers, V. Rosenhaus, and M. Smolkin. in progress. 

[80] F. Larsen and F. Wilczek, “Renormalization of black hole entropy and of the gravitational 
coupling constant,” Nucl.Phys. B458 (1996) 249-266, arXiv:hep-th/9506066 [hep-th]. 

[81] H. Osborn and G. Shore, “Correlation functions of the energy momentum tensor on spaces 
of constant curvature,” Nucl.Phys. B571 (2000) 287-357, arXiv :hep-th/9909043 
[hep-th] . 

[82] L.-Y. Hung, R. C. Myers, and M. Smolkin, “Some Calculable Contributions to Holographic 
Entanglement Entropy,” JHEP 1108 (2011) 039, arXiv: 1105.6055 [hep-th]. 

[83] D. Marolf, D. Minic, and S. F. Ross, “Notes on space-time thermodynamics and the 
observer dependence of entropy,” Phys.Rev. D69 (2004) 064006, arXiv: hep-th/0310022 
[hep-th] . 

[84] J. Bhattacharya, M. Nozaki, T. Takayanagi, and T. Ugajin, “Thermodynamical Property of 
Entanglement Entropy for Excited States,” Phys.Rev.Lett. 110 no. 9, (2013) 091602, 
arXiv:1212.1164. 

[85] D. D. Blanco, H. Casini, L.-Y. Hung, and R. C. Myers, “Relative Entropy and 
Holography,” JHEP 1308 (2013) 060, arXiv: 1305.3182 [hep-th], 

[86] G. Wong, I. Klich, L. A. Pando Zayas, and D. Vaman, “Entanglement Temperature and 
Entanglement Entropy of Excited States,” JHEP 1312 (2013) 020, arXiv: 1305.3291 
[hep-th] . 

[87] D. Allahbakhshi, M. Alishahiha, and A. Naseli, “Entanglement Thermodynamics,” JHEP 
1308 (2013) 102, arXiv:1305.2728 [hep-th], 

[88] M. Smolkin and S. N. Solodukhin, “Correlation functions on conical defects,” Phys.Rev. 
D91 no. 4, (2015) 044008, arXiv: 1406.2512 [hep-th]. 

[89] H. Osborn and A. Petkou, “Implications of conformal invariance in field theories for general 
dimensions,” Annals Phys. 231 (1994) 311-362, arXiv:hep-th/9307010 [hep-th]. 

[90] E. Bianchi and R. C. Myers, “On the Architecture of Spacetime Geometry,” 

Class.Quant.Grav. 31 no. 21, (2014) 214002, arXiv: 1212.5183 [hep-th]. 

[91] S. Forte and J. I. Latorre, “A Proof of the irreversibility of renormalization group flows in 
four-dimensions,” Nucl.Phys. B535 (1998) 709-728, arXiv:hep-th/9805015 [hep-th]. 

[92] L. S. Brown and J. C. Collins, “Dimensional Renormalization of Scalar Field Theory in 
Curved Space-time,” Annals Phys. 130 (1980) 215. 

[93] L. S. Brown, “Stress Tensor Trace Anomaly in a Gravitational Metric: Scalar Fields,” 
Phys.Rev. D15 (1977) 1469. 


- 37 - 


[94] R. Camporesi, “The Spinor heat kernel in maximally symmetric spaces,” 
Commun.Math.Phys. 148 (1992) 283-308. 

[95] A. Higuchi, “Symmetric Tensor Spherical Harmonics on the N Sphere and Their 
Application to the De Sitter Group SO(AT,l),” J.Math.Phys. 28 (1987) 1553. 

[96] A. Cappelli, D. Friedan, and J. I. Latorre, “C theorem and spectral representation,” 
Nucl.Phys. B352 (1991) 616-670. 

[97] A. C. Wall, “Ten Proofs of the Generalized Second Law,” JHEP 0906 (2009) 021, 
arXiv:0901.3865 [gr-qc] . 

[98] D. V. Fursaev and S. N. Solodukhin, “On the description of the Riemannian geometry in 
the presence of conical defects,” Phys.Rev. D52 (1995) 2133-2143, arXiv:hep-th/9501127 
[hep-th]. 

[99] P. Calabrese and J. L. Cardy, “Entanglement entropy and quantum field theory: A 
Non-technical introduction,” Int.J. Quant.Inf. 4 (2006) 429, arXiv: quant-ph/0505193 
[quant-ph] . 

[100] D. Vassilevicli, “Heat kernel expansion: User’s manual,” Phys.Rept. 388 (2003) 279-360, 
arXiv:hep-th/0306138 [hep-th]. 

[101] D. V. Fursaev, A. Patrushev, and S. N. Solodukhin, “Distributional Geometry of Squashed 
Cones,” Phys.Rev. D88 no. 4, (2013) 044054, arXiv: 1306.4000 [hep-th]. 

[102] A. Bhattacharyya, M. Sharma, and A. Sinha, “On generalized gravitational entropy, 
squashed cones and holography,” JHEP 1401 (2014) 021, arXiv: 1308.5748 [hep-th], 

[103] X. Dong, “Holographic Entanglement Entropy for General Higher Derivative Gravity,” 
JHEP 1401 (2014) 044, arXiv: 1310.5713 [hep-th], 

[104] J. Camps, “Generalized entropy and higher derivative Gravity,” JHEP 1403 (2014) 070, 
arXiv:1310.6659 [hep-th]. 

[105] A. Bhattacharyya and M. Sharma, “On entanglement entropy functionals in higher 
derivative gravity theories,” JHEP 1410 (2014) 130, arXiv: 1405.3511 [hep-th]. 

[106] J. Dowker, “Quantum Field Theory on a Cone,” J.Phys. A10 (1977) 115-124. 

[107] J. Dowker, “Vacuum Averages for Arbitrary Spin Around a Cosmic String,” Phys.Rev. D36 
(1987) 3742. 

[108] S. Deser and R. Jackiw, “Classical and Quantum Scattering on a Cone,” 
Commun.Math.Phys. 118 (1988) 495. 


- 38 - 


